Localization and critical diffusion of quantum dipoles in two dimensions. 
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We discuss quantum propagation of dipole excitations in two dimensions. This problem differs 
from the conventional Anderson localization due to existence of long range hops. We found that the 
critical wavefunctions of the dipoles always exist which manifest themselves by a scale independent 
diffusion constant. If the system is T-invariant the states are critical for all values of the parameters. 
Otherwise, there can be a "metal-insulator" transition between this "ordinary" diffusion and the 
Levy- flights (the diffusion constant logarithmically increasing with the scale). These results follow 
from the two-loop analysis of the modified non-linear supermatrix cr-model. 

PACS numbers: 71.23.-k, 71.55.Jv 
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Anderson [1] showed that a quenched disorder can lo- 
cahze a quantum particle, i.e. completely suppress its 
diffusion. Later [2, 3], it was realized that in two di- 
mensions {2D) localization occurs for an arbitrary weak 
disorder. This conclusion was reached by studying the 
scaling behavior of the dimensionless Thouless conduc- 
tance g{L) as a function of the linear size of the system 
L (the observable electrical conductance of the system 
of the e charged particles is given by 5 x e^/S, and we 
will set Planck constant h = 1 hereinafter) . Localization 
implies that g/, — >■ as L — > 00. This is always true when 
the time reflection symmetry (T-invariance) is broken (so 
called unitary ensemble, GUE). For T-invariant systems 
it is still correct if the orbital and spin degrees of freedom 
are decoupled or when the particles have an integer spin 
(orthogonal ensemble, GOE). For the particles with half- 
integer spin, the theory [4] predicts that the spin-orbital 
coupling causes antilocalization g{L — >■ 00) — cxd if disor- 
der is weak, while for a stronger disorder g{L — > 00) — >■ 
(metal-insulator transition for symplectic ensemble) [5] . 

Besides the current carrying charged particles impor- 
tant objects in the many-body theory are neutral exci- 
tations (NEX) - the bound states of two particles with 
opposite charges. One can name excitons in semiconduc- 
tors, optical phonons in polar crystals, dipole excitations 
in granular superconductors, vacancy-interstitial excita- 
tions in Wigner crystals or vortex lattices, etc. Being the 
lowest energy excitations, the NEX not only determine 
the low temperature energy transport but also provide a 
thermal bath for more energetic charge excitations, thus 
influencing the charge transport. As for any neutral par- 
ticles, the number of NEX is not conserved {e.g. electron 
and hole can annihilate each other) . Each NEX has a fi- 
nite energy and cannot simply disappear. However, their 
number non-conservation facilitates long range hops me- 
diated by virtual photons (this leads also to dipole-dipole 
interactions). In this Letter we investigate the effect of 
the long-range hops on the localization of NEX. 

We will be interested in 2D quantum dipoles - NEX 



whose annihilation [creation] operators ^^(r) [&J,(r)] are 
characterized by the additional index a = x, z. The 
pair {bx,bz) transforms under rotations similar to a vec- 
tor in 2D plane {x,z). For a small density of NEX, 
we can neglect the interaction between them and use 
a bilinear form of ha , as the Hamiltonian T-L = 
/ d?Vi(Pv2b]^{vi)Hafj{Yi - Y2)bfj{v2), (wc imply summa- 
tion over the repeated indices a,j3 = x,z hereinafter). 
The absorption-emission of the virtual photons (with in- 
finite speed) results in a long-range hopping term [6] 
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A(5„^|rp-2r,r^)/(27r|r|4), 
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additional to the local Hamiltonian 



H: 
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J(ri-r2)[i?°^(-*V) + K;5(r)]. (2a) 



Constant A in Eq. (1) encodes dipole-photon transitions 
matrix elements and the energy denominators. 

The hops (1) qualitatively modify the Anderson lo- 
calization. We derived and analyzed the renormalization 
group (RG) equations describing the scaling of the Thou- 
less conductance g{L, A). For the GOE we found the sta- 
ble fixed line 5c (A) on the g — X plane, which corresponds 
to the the critical state, see Fig. la. Remarkably, (/c(A) 
may be large, which makes this line accessible in the per- 
turbativc RG. For the GUE, we discover the unstable line 
5mit(A) separating the antilocalization and the localiza- 
tion behaviors, see Fig. lb. We argue that the latter is 
the precursor of the critical state similar to the one in the 
GOE, which, unfortunately, occurs at (7 < 1 and thus, is 
out of reach of the perturbative RG analysis. 

The quantum motion of the dipoles in the rotationally 
symmetric clean system is described by 



2mi(k) ' 2m2(k) 



(2b) 



where mi_2 (k) are analytic functions of the momentum 
k, invariant with respect to the lattice symmetry group. 
The Pauli matrices a^, i = x,y,z act in 2 x 2 space of 



FIG. 1. RG flows for (a) GOE and (b) GUE, and A{X) < 1 
is defined in Eq. (9). Sofid lines are accessible by the pertur- 
bative RG, and dashed ones are conjectures. 



dipole components xz. The last term m Eq. (2b) breaks 
the T-invariance by the magnetic field oc h (The Hamil- 
tonian of NEX cannot contain the vector potential) . At 
h = 0, we can diagonalizc ff" + iJ''' (see Fig. 2a): 



2mi{k)' 



i?+(k)=i?_(k)- 



2m2(k) 



+ A. (2c) 



Note that E± (k) are analytic functions of k, even though 
H''^ lifts the degeneracy at fc = protected by symmetry. 

The 2x2 matrix V{r) in xz-space is a Gaussian dis- 
order breaking all the system spatial symmetries, 



f (r) ® V{r')) = Sir - r') ^ 



(2d) 



i—0,x.y,z 



where (Tq = 1. For T-invariant systems, Uy = 0. The 
rotational symmetry after the disorder averaging requires 
Ux = Uz- The elastic mean free time r is given by 1/t = 
TTi^^j^Q ^ u|, with the total density of states being 



^W=E / '5[e-i?±(k)]dV(2^)'- (3) 



We assume that the disorder is weak, er ^ 1. 

For the specific form of Eq. (1), the time evolution 

U{t) = exp(it / <fvi(Pr2bi{vi)H'^p{Y^2)h{^2)) can be 
described locally by introducing additional fields 



U^''{t) = j 2?a(r)Pa*(r)exp (^-^ J cPrcHdi 



1 



(4) 
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Due to the locality of the Poisson equation describing 
virtual photons, operator M is differential rather than 
integral one and this enables us to develop a renormaliz- 
able theory of localization. 

To understand the effect of the long hops (1) on the 
localization, consider the dipole with high energy, e A, 
see Fig. 2a. For the clean system, the wave functions are 
comprised by plane waves with wavevectors k±(t) and 
oscillate rapidly. The long-range hops are thus irrelevant. 





FIG. 2. a) Band structure for the clean system; b) Long hop 
assisted by the impurities (crosses); c) Semiclassical paths for 
one-loop weak localization; d) Long range hops and semiclas- 
sical paths for two-loop renormalization; e) Two loop semi- 
classical paths including only local evolution. 



In the presence of disorder a dipole is characterized by 
the probability density n(r) of finding it at a point r. 
Under the assumption of Markovian time evolution. 



dtn{r) = j <f-RW{K) [n[r + R) - n(r)] 



(5a) 



At large enough distances and times, Eq. (5a) reduces to 
the diffusion equation with diffusion constant D deter- 
mined by the rates kernel VF(R): 



(5b) 



To estimate M^(R), consider the transition between 
wave-packets 0q , formed by plane- waves close to the 
energy shell and normalized / d^r^a'' [^q"*]* = 1, see 
Fig. 2b. The disorder can e.g. scatter 0a •* into a vir- 
tual small momentum state, this state can be transferred 
far by the Hamiltonian (1), and then scattered by disor- 
der back onto the energy shell (upper part of Fig. 2b). 
The amplitude for this process, J(R), is given by 

A (Rii.i"(.f ')• - 2t,i'lR„i!j(»f )•) 

2im' • 

where wi^'^ = / d^r,/)^^' V"'". We substitute Eq. (6) into 
the golden rule formula, average it over the disorder (2d) , 
and use the normalization of (b'£' . We find 



I^,,(R) = 27ri.(| J(R)p) = AV (47r3i.eV2|i?|4) , (7) 
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Substituting Eq. (7) into Eq. (5b) and taking into account 
only hops with L < \R\ < L + 6L, we obtain a positive 
correction to the Thouless conductance g = vD: 



= <5M/(27r2); ^ = ln(L/io), 



(8) 



where Lq is a microscopic scale, and [see Eq. (7)] A = 
A^/ (4r^e'') ; re » 1, e » A. For e in resonance with E±, 
the summation of all orders of perturbation theory in A 
gives the Breit-Wigner type formula 
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2^2 



[4 + T'-e+] [4 + T^e] 



<1, e± = e-i?±(0), (9) 



i.e. even at A > 1 (e ^ A), the contribution of the long- 
range hops cannot exceed the unitary limit. 

Levy flight term (8) is not the only logaritmical con- 
tribution to the conductance. The interference of close 
time-reversed paths {h,Uy = 0) controlled by local term 
}{sr (Yig^ 2c) yields the weak localization correction [3] 



Sg^i = -Se/ (2^2) 



(10) 



Since A < 1, long hops (8) never overcome weak local- 
ization but can almost compensate it if \A — 1| ^ 1. 

In this case, it is worthwhile to evaluate the two-loop 
contribution. When included as intermediate steps into 
interference contribution of Fig. 2c, the long range hops 
simply change the bare diffusion constant D. However, 
Sg^i does not contain D at all, i.e. irreducible inter- 
ference processes are important. One of them. Fig. 2e, 
vanishes as the path (1) interferes destructively with path 
(2) and constructively with paths (3), (4). It is the long 
hop part of the Hamiltonian, Fig. 2d, that leads to the 
logarithmic correction to g. Modified non-linear cr-model 
described below yields a new RG equation for GOE 



dg 
di 



A~l 



STT^g 



+ 0(l/.g2), (11) 



yielding the RG flow of Fig. la and the stable fixed line 



g, (A) = A/ [4^2 (1 _ ^)] 



(12) 



Note, that gc ^ 1 in the limit 1 — ^ ^ 1 i.e. the two 
loop approximation is sufficient. At 5 < 1, Eq. (12) is not 
applicable. However, there is a strong reason to believe 
that the critical line terminates at g = 0, A = point. 
Indeed, ii g < 1, the local part of the evolution can pro- 
duce only small contributions Sg ~ S£glng. Arguments 
leading to the estimate (8) remain valid. Indeed, one 
can chose the exact localized eigenfunctions of the Hsr 
instead of the wave packets and obtain 



dg/d£ = A/(27r2) + 0{glng) > 0, .9 < 1, 



(13) 



i.e., g = is unstable for any A > 0. As dg/d£ < for 
g ^ 1, see Eq. (11), the fixed line has the form of Fig. la. 



Application of the magnetic field removes the time re- 
versal symmetry, h,Uy 7^ 0, and thus suppresses the con- 
tribution of time reversal paths of Fig. 2 c,d at distances 
larger than Lh = ^/D/uJh, where cOh — h'^T + ttvu^. The 
interference between path (1) and paths (3,4) is also de- 
stroyed and there is no more cancelation of the local in- 
terferences in the second loop. We find for GUE 



dg_ 

di 



A 

2^2 



(14) 



similarly to results of Ref. [7] for an electron in random 
magnetic field with long-range correlations. 
Equation (14) has an unstable fixed line 



g™t (A) = {in^Ay 



(15) 



corresponding to the metal-insulator transition. This re- 
sult is controllable for A < 1. The RG flow of Fig. 2b 
follows from Eq. (15) and Eq. (13) for g < I. 

Let us discuss the relation of our results to earlier 
works. At first glance, the long jumps are equivalent 
to artificial models with long-range random links known 
as random band matrices, REM [8, 9], (the matrix ele- 
ments are the Gaussian variable with (hfj) ~ l/|i ~ 
{hij) — 0) The dipole hop model (1) is different as it 
contains an infinite number of the long hop loops {e.g. 
the product of the matrix elements (/i 12 ^23/134/141) for 
non-coinciding points 1,...,4 vanishes for REM and it 
is finite for the dipolar interaction) and determines all 
of the interference contributions of Fig. 2d. The quan- 
tum dipole problem with long range hops was considered 
by Lcvitov [10]. He started from the strongly localized 
states, g — > 0, and replaced the strongly localized pairs of 
Fig. 2b by one effective dipole on each linear scale. This 
approach misses the logarithmic contributions of the mul- 
tiple short range hops. We believe that in such a way one 
can show the instability of g = state but can not ob- 
tain the critical line (g oc A claimed in Ref. [11] is neither 
Thouless conductance nor ac-conductancc. Both those 
quantities are oc A^, see Eq. (13)). 

Supersymmetric non-linear a (nla) model is usually 
formulated [12] in terms of the 8-component supervec- 
tor. The 8-dimensional space where this vector resides 
can be presented as a direct product of 3 two-dimensional 
subspaces, RA, N and 5, for the retarded-advanced, 
Gorkov-Nambu, and the fermion-boson sectors. The ma- 
trix structure of the Hamiltonians (1) and (2), brings up 
an additional two dimensional subspace (without sub- 
script), so the resulting supervector b is 16-component. 
One more supervector a with the same structure as b is 
introduced to decouple the long-range hops as it is shown 
in Eq. (4). Those two vectors can be united in one 32- 
dimensional supervector ip'^ = (a^ ,b'^)°'^ . The disorder 
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averaging is then performed as 



•RA 



RA 



N 



) a 

(a- 



) Iff - 



(16a) 



(...)= / ••■exp 

k = k ® t"-" ® t] A 
C = (7(8)1"'' (8)1; (7=1 

where ... in the LHS stands for any combina- 
tion of advanced/retarded Green functions G^'^ = 

eTi^l'^- Ho-V - HirT io) ^ and • • • in the RHS 
for the corresponding sources [12], whose form is not im- 
portant for RG. Term Lq in the Lagrangian L = Lo+Lint 
describes the property of the clean system 



110 = 



iXM/2- 



AO^ 



ibIHIo(-iV)b| 









'k^ ® k" 






(8 1 + 


\2mi 




2m2 



ha, 



V 



-V„(7„ 1 



where / EE 1^^®1^^®19, and S3 = 1^^ ® o-f (8 1^*. The 



term Lint originates from the disorder averaging (2d): 

iha^bfSY. (16c) 



Lint W\^\ ' 



As for the conventional nZcr-modcl the interaction is 
decoupled by introduction the additional field Q (xh®h 
and fixing = / in the saddle point approximation is 
valid as long as er 3> 1. For a = 0, the field b for a fixed 
Q is massive and can be integrated out producing con- 
vergent gradient expansion in nla model. However, field 
a is massless and its fluctuations lead to the logarithmic 
corrections discussed above. Therefore, only "fast" field 
■0 with the momenta close to the energy shell k±{e) can 
be included into the effective theory for Q, whereas fields 
with the smaller momenta must be kept in the theory: 



(•••) = 



(17a) 



Q = CQ^C^; Q 



it _ 
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RA 



1 



N 



Here, ip is the 32-componcnt smooth on the scale of 
the mean free path field satisfying the constraints of 
Eq. (16b); Q is the smooth 8x8 supermatrix [12]. The 
entries in Eq. (17a) are 



F = -str 



2iukQ — LOh ( S3 (5 



C 



y 



d2rV'(r)B(g)V'(r) 



(17b) 



el 



2r 



1 + ha„ 



where the bare conductance is given by 
d^k TdE±dE± 
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± 



(27r)' 2 dk^ dk^ 



<5[6-i?±(k)]. 



Note that, in contrast to non-local RBM models [8], 
the theory (17) is renormalizable. Indeed, F[Q] of 
Eq. (17a) includes only relevant terms allowed by = I. 
The same applies to the Lagrangian C as the natural di- 
mensions of fields [Va] = [b] = 1 while a itself can not 
enter: C is invariant with respect to constant shift of a. 

We checked that the fluctuations of (^-matrix do not 
change the Lagrangian C and, thus, the coefficient A in 
Eqs. (11)-(14). This relates to the fact that the averaged 
density of states cannot have corrections from mixing of 
retarded and advanced sectors described by Q matrix. 

In order to sum up leading logarithmic diver- 
gences, we performed two-loop RG analysis of theory 
(17b), see Ref. [13] for details. We represent Q = 



(16b) V 



{l + iP)/{l~iP) 



V and -0 



Vipy , where 



'0<, V are slow variables and the other are fast. We inte- 
grated out the fast variables in the second loop approxi- 
mation to obtain the effective free energy for slow ones. 
Varying the result of the integration with respect to a 
gauge invariant cutoff we obtain Eqs. (11) and (14). 

In conclusion, we performed the scaling analysis of the 
localization problem of 2D quantum dipoles, see Fig. 1. 
For T-invariant systems, the Thouless conductance tends 
to a finite value. Breaking the T-invariance leads to the 
transition between the antilocalization and a critical be- 
havior. To describe those phenomena we constructed a 
novel version of the non-linear a-model, Eq. (17), and 
derived the second loop RG equations. This model also 
allows studies [13] of the multifractal properties of the 
wave- function using methods of Ref. [12]. 
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